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^^ , ABSTRACT: We study the non stationary solutions of Fokker-Planck equations associated to 

Mh, either stationary or non stationary quantum states. In particular we discuss the stationary states 

f—t . of quantum systems with singular velocity fields. We introduce a technique that allows arbitrary 

2 , evolutions ruled by these equations to account for controlled quantum transitions. As a first signfi- 

^^, cant application we present a detailed treatment of the transition probabilities and of the controlling 

^ , time-dependent potentials associated to the transitions between the .stationary, the coherent, and 

K^ , the squeezed states of the harmonic oscillator. 

1. Introduction 

In a few recent papers [1] the analogy between diffusive classical systems and quantum systems 
has been reconsidered from the standpoint of the stochastic simulation of quantum mechanics [2] , 
[3] , [4] and particular attention has been devoted there to the evolution of the classical systems 
associated to a quantum wave function when the conditions imposed by the stochastic variational 
principle are not satisfied (non extremal processes). The problem studied in those papers was 
the convergence of an arbitrary evolving probability distribution, solution of the Fokker-Planck 
equation, toward a suitable quantum distribution. In [1] it was pointed out that, while the correct 
convergence is achieved for a few quantum examples, these results cannot be considered general 



as shown in sonic counterexamples: in fact not only for particular non stationary wave functions 
(as for a minimal uncertainty packet), but also for stationary states with nodes one does not 
recover in a straightforward way the correct quantum asymptotic behaviour. For stationary states 
with nodes the problem is that the corresponding velocity field to consider in the Fokker-Planck 
equation shows singularities at the locations of the nodes of the wave function. These singularities 
effectively separate the available interval of the configurational variables into non communicating 
sectors which trap any amount of probability initially attributed and make the system non ergodic. 

In a more recent paper [5] it has been shown that for transitive systems with stationary 
velocity fields (as, for example, a stationary state without nodes) we always have an exponential 
convergence to the correct quantum probability distribution associated to the extremal process, 
even if we initially start from an arbitrary non extremal process. These results can also be extended 
to an arbitrary stationary state if we separately consider the process as confined in every region of 
the configuration space between two subsequent nodes. 

In the same paper [5] it has been further remarked that while the non extremal processes should 
be considered virtual, as the non extremal trajectories of classical Lagrangian mechanics, they can 
however become physical, real solutions if we suitably modify the potential in the Schrodinger 
equation. The interest of this remark lies not so much in the fact that non extremal processes are 
exactly what is lacking in quantum mechanics in order to interpret it as a totally classical theory 
of stochastic processes (for example in order to have a classical picture of a double slit experiment 
[6] ) , but rather in the much more interesting possibility of engineering and controlling physically 
realizable evolutions of quantum states. This observation would be of great relevance, for instance 
to the study and the description of (a) transitions between quantum states (b) possible models for 
quantum measurements [3] and (c) control of the dynamics of quantum-like systems (for instance 
charged beams in particle accelerators) [7]. 

In particular, case (c) is being studied in the framework of Nelson stochastic mechanics which 
is an independent and self-consistent reformulation of quantum mechanics [2], [3] and can be 
applied in other areas of physical phenomenology. For instance it can usefully account for systems 
not completely described by the quantum formalism, but whose evolution is however strongly 
influenced by quantum fluctuations, i.e. the so-called mesoscopic or quantum-like systems. This 
behaviour characterizes, for example, the beam dynamics in particle accelerators and there is 
evidence that it can be described by the stochastic formalism of Nelson diffusions [1], [7] since 
in these quantum-like systems, trajectories and transition probabilities acquire a clear physical 
meaning, at variance with the case of quantum mechanics. 

On the other hand, quantum behaviours can be simulated by means of classical stochastic 
processes in a by now well defined and established framework [2]. A stochastic variational prin- 
ciple provides a foundation for that, in close analogy with classical mechanics and field theory 
[3]. In this scheme the deterministic trajectories of classical mechanics are replaced by the random 
trajectories of diffusion processes in configuration space. The programming equations derived from 
the stochastic variational principle are formally identical to the equations of the Madelung fluid [8] , 



the hydrodynamical equivalent of Schrodinger equation [9]. On this basis, it is possible to develop 
a model whose phenoincnological predictions coincide with those of quantum mechanics for all 
the experimentally measurable quantities. Within this interpretative code stochastic mechanics is 
nothing but a quantization procedure, different from the canonical one only formally, but com- 
pletely equivalent from the point of view of the physical consequences: a probabilistic simulation of 
quantum mechanics, providing a bridge between this fundamental physical theory and stochastic 
differential calculus. However, it is well known that the central objects in the theory of classical 
stochastic processes, namely the transition probability densities, seldom play any observable role in 
stochastic mechanics and must be considered as a sort of gauge variables. Several generalizations 
of Nelson stochastic quantization have been recently proposed to allow for the observability of 
the transition probabilities: for instance, stochastic mechanics could be modified by means of non 
constant diffusion coefficients [1]; alternatively, it has been suggested that the stochastic evolution 
might be modified during the measurement process [10]. 

The aim of the present paper is instead to show how the transition probabilities associated to 
Nelson diffusion processes can play a very useful role in standard quantum mechanics, in particular 
with regard to describing and engineering the dynamics of suitably controlled quantum evolutions 
and transitions. More precisely, we consider the following problem in the theory of quantum control: 
given an initial probability distribution pi associated to an arbitrarily assigned quantum state ipi, 
we study its time evolution with the drift associated to another arbitrarily assigned quantum state 
Tpf, to determine the controlling time-dependent potential Vc{x,t) such that, I) at any instant of 
time the evolving probability distribution is that associated to the wave function solution of the 
Schrodinger equation in the potential Vc{x,t), and that, II) asymptotically in time the evolving 
distribution converges to the distribution pf associated to ipf. 

After introducing the formalism of Nelson stochastic mechanics to describe quantum evolutions 
in Sections 2 and 3, we provide in Sections 3 and 4 a self-contained review of the Sturm-Liouville 
problem for the Fokker-Planck equation and the techniques of solution for the Nelson diffusions 
associated both to nonstationary and stationary quantum states. In Section 5 we discuss in detail 
the example of the harmonic oscillator, explicitely solving for the transition probability densities 
of the ground and of the low lying excited states. Sections 6, 7 and 8 are devoted to the study and 
the solution of the problem outlined above, discussing the potentials associated to the definition 
of controlled quantum evolution, and modelling transitions. Two explicit examples are studied in 
detail: the controlled transition between the invariant probability densities associated to the ground 
and the first excited state of the harmonic oscillator, and the controlled evolution between pairs of 
coherent or squeezed wave packets. In these cases the problem can be solved completely, yielding 
the explicit analytic form of the evolving transition probabilities and of the evolving controlling 
potentials at all times. Finally, in Section 9 we present our conclusions and discuss possible future 
extensions and applications of the technique introduced in the present paper, with regard to the 
discussion of anharmonic quantum and quantum-like systems, the role of instabilities in the initial 
conditions, and the implementation of optimization procedures. 



2. Fokker-Planck equations and quantum systems 

Here we will recall a few notions of stochastic mechanics in order to fix the notation. The con- 
figuration of a classical particle is promoted to a vector Markov process ^(t) taking values in R"^. 
This process is characterized by a probability density p{r, t) and a transition probability density 
p(r,t| r',t') and its components satisfy an Ito stochastic differential equation of the form 

^OW = t'(+)j(eW,i)rfi + dfi.j{t) , (2.1) 

where i'(+)j are the components of the forward velocity field. However here the fields f (+)j are 
not given a priori, but play the role of dynamical variables and are consequently determined by 
imposing a specific dynamics. The noise ri(t) is a standard Wiener process independent of ^(t) and 
such that 

Et{drjj{t)) ^ , I]t{dr]j{t)dr]k{t)) ^2D5jkdt, (2.2) 

where drj(t) = rj(t + dt) —i]{t) (for dt > 0), D is the diffusion coefficient, and Et are the conditional 
expectations with respect to £,{t). In what follows, for sake of notational simplicity, we will limit 
ourselves to the case of one dimensional trajectories, but the results that will be obtained can be 
immediately generalized to any number of dimensions. We will suppose for the time being that 
the forces will be defined by means of purely configurational potentials, possibly time-dependent 
V{x,t). A suitable definition of the Lagrangian and of the stochastic action functional for the 
system described by the dynamical variables p and t;(+) allows to select, the processes which 
reproduce the correct quantum dynamics [2], [3]. In fact, while the probability density p{x,t) 
satisfies, as usual, the forward Fokker-Planck equation associated to the stochastic differential 
equation (2.1) 

dtp = Ddlp - d^{v(+)p) ^ d^{Dd^p - V(+)p) , (2.3) 

the following choice for the Lagrangian field 

L{x, t) = y w^+) (x, t) + mDa,«(+) {x, t) - V{x, t) , (2.4) 

enables to define a stochastic action functional 

/■ti 
A= / E[L(C(t),i)]rfi, (2.5) 

which leads, through the stationarity condition 5A = 0, to the equation 

a,5+IMl! + F-2mi?2^VP^0 (2.6) 

2m y/p 

The field S{x,t) is defined as 

S{x, t) = -J\ [L{as), s) I m =x]ds + E [Si (e(ti)) I m = x] , (2.7) 
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where S'i(-) = S{-,ti) is an arbitrary final condition. By introducing the function R{x,t) = 
y/p{x,t) and the dc Broglic Ansatz 

i;{x, t) = i?(x, t) e'Si^^t)/2rnD ^ ^2.8) 

equation (2.6) takes the form 

dtS+^-^^ + V -2mD'^^^0, (2.9) 

and the complex function ip satisfies the Schrodinger-like equation 

i{2mD)dttp = Hip = -2mD'^dltp + V^p. (2.10) 

If the diffusion coefficient is chosen to be 

D=i^, 2.11 

Zm 

we recover exactly the Schrodinger equation of quantum mechanics. Different choices of D allow 
instead to describe the effective quantum-like dynamics of more general systems. 

On the other hand, if we start from the (one-dimensional) Schrodinger equation (2.10) with the 
de Broglie Ansatz (2.8) and the diffusion coefficient (2.11), separating the real and the imaginary 
parts as usual in the hydrodynamical formulation [8], we recover equations (2.3) and (2.6) with 
p ^ R^ = j-f/ip and the forward velocity field 

V(^+^{x,t)^-d,S+^d^(\nR^). (2.12) 

m Im 



3. The Sturm Liouville problem and the solutions of the Fokker Planck equation 

Let us recall here (see for example [11]) a few generalities about the techniques of solution of 
the Fokker-Planck equation with D and vi^\ two time-independent continuous and differentiable 
functions defined for x G [a, h] and t > to, such that D{x) > 0, and W(+)(x) has no singularities in 
(a, 6). The Fokker-Planck equation then reads 

dtp = dliDp) - a,(«(+)p) = d, [d,{Dp) - «(+)p] . (3.1) 

The conditions imposed on the probabilistic solutions are of course 

p{x, t) > , a<x<b,tQ<t, 

rb (3.2) 

p{x, t) dx ^ I , to < t , 

and from the form of (3.1) the second condition takes the form 

[d,{Dp)-V(^+^p]^,^ = 0, to<t. (3.3) 
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Suitable initial conditions will be added to produce the required evolution: for example the tran- 
sition probability density p{x, i|xo, to) will be selected by the initial condition 



lim p{x,t) = pt„{x) = 5{x ~ xo) . 



(3.4) 



It is also possible to show by direct calculation that 



h(x) = 7V"^e"/'"°'^'^''"''<+>^'^'''^'°'''^'''''^ 



N= / e 



^- J[D'{x)-Vf^^{x)]/D(x)dx j^^ \ /g gx 



is always an invariant (time independent) solution of (3.1) satisfying the conditions (3.2) (here the 
prime symbol denotes differentiation). One should observe however that relation (3.1) is not in 
the standard self-adjoint form [12]; this fact notwithstanding, if we define the new function g{x, t) 
by means of 

p{x,t)^ ^/h(x)g{x,t), (3.6) 

it is easy to show that g{x,t) obeys an equation of the form 

dtg = Cg , 



where the operator C acting on positive normalizable functions ip{x) and defined by 



C(p 



dx 



r{x) 



d(p{x) 



dx 



q{x)ip{x) , 



with 



r{x) =D{x) > 0, 

^(^)^i^'(-) -«(+)(-)]' [D'{x)-v,^,{x)]' 



(3.7) 



(3.8) 



(3.9) 



AD{x) 2 

is now self-adjoint. By separating the variables by means of g{x, t) — '-f{t)G{x) we have j{t) — e^^* 
while G must be solution of a typical Sturm-Liouville problem associated to the equation 



CG{x)+XG{x) =0, 



(3.10) 



(3.11) 



with the boundary conditions 

[D'{a) - 'y(+)(a)] Gia) + 2D{a)G'{a) = , 
[D'{b) - i;(+)(6)]G(6) + 2D{h)G'{b) = . 

It is easy to see that A = is always an eigenvalue for the problem (3.10) with (3.11), and that 
the corresponding eigenfunction is ^Jh{x) as defined from (3.5). 

For the differential problem (3.10) with (3.11) we have that [12] the simple eigenvalues A„ 
will constitute an infinite, monotonically increasing sequence and the corresponding eigenfunction 
Gn{x) will have n simple zeros in (a, b). This means that Ao — 0, corresponding to the eigenfunction 



Gq{x) — y^h(x) which never vanishes in (a, 6), is the lowest eigenvalue, and that all the other 



eigenvalues are strictly positive. Moreover the eigenfunctions will form a complete ortlionormal set 
of functions in L'^([a,b]] [13]. As a consequence, the general solution of equation (3.1) satisfying 
the conditions (3.2) will have the form 

oo 

p(x,t)-^c„e-^"*yM^G„(x), (3.12) 

n=0 

with Co = 1 for normalization (remember that Aq = 0). The coefficients c„ for a particular solution 
are selected by an initial condition 

p(x,tj) = Pj„(x), (3.13) 

and are then calculated from the orthonormality relations as 

c„ = / pin{x) —^== dx . (3.14) 

Ja \/h{x) 

In particular for the transition probability density we have from (3.4) that 

G„(xo) 



\/h(xQ) 



(3.15) 



Since Ao = and A„ > for n > 1, the general solution (3.12) of (3.1) has a precise time evolution: 
all the exponential factors in (3.12) vanish as t ^ +oo with the only exception of the term n = 
which is constant, so that exponentially fast we will always have 



lim p{x,t) = coy^h{x)Go{x) ~ h{x) . (3.16) 

t — ^ + oo 

Therefore the general solution will always relax in time toward the invariant solution h{x). As a 
consequence the eigenvalues A„ which solve the Sturm Liouville problem define the physical time 
scales of the decay. By the structure of equations (3.7)-(3.11) we see that tuning the choice of 
the physical parameters that enter in the diffusion coefficient and in the forward velocity field 
allows for different sets of eigenvalues which define different sets of time scales. Hence, the rate of 
convergence can be fixed as to yield fast decay, slow decay, or even, on proper observational scales, 
quasi metastable behaviours, according to what kind of physical evolution between quantum states 
one wants to realize. This point will be further discussed and elucidated in Section 6. 



4. Processes associated to stationary quantum states 

Let us consider now a Schrodinger equation (2.10) with a time-independent potential V{x) which 
gives rise to a purely discrete spectrum and bound, normalizable states. Let us introduce the 
following notations for stationary states, eigenvalues and eigenfunctions: 

M^,t)^cj,n{x)e-'''-"'', 

n" „ (4-1) 
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Taking into account relation (2.11), the previous eigenvalue equation can also be recast in the 
following form 

i?C--^^^0«. (4.2) 

For these stationary states the probability densities are the time-independent, real functions 

Pn{x)^\lbn{x,t)\^^4'l{^), (4.3) 

while the phase and the amplitude of '0„ are from (2.8) 

Snix,t)^-Ent, Rn{x,t)^(j)n{x), (4.4) 

SO that the associated velocity fields are from (2.12) 

Each f (+)„ is time-independent and presents singularities in the nodes of the associated eigenfunc- 
tion. Since the n-th eigenfunction of a quantum system with bound states has exactly n simple 
nodes xi, . . . ,a;„, the coefficients of the Fokker-Planck equation (2.3) are not defined in these n 
points and it is necessary to solve it in separate intervals by imposing the correct boundary con- 
ditions connecting the different sectors. In fact these singularities effectively separate the real axis 
in n + 1 sub-intervals with walls impenetrable to the probability current. Hence the process will 
not have an unique invariant measure and will never cross the boundaries fixed by the singularities 
of W(+)(x): if the process starts in one of the sub-intervals, it will always remain there [14]. 

As a consequence, the normalization integral (3.2) (with a = — oo and b — +oo) is the sum of 
n + 1 integrals over the sub-intervals [xk, Xk+i] with fc = 0, 1, . . . , n (where we understand, to unify 
the notation, that xq = — oo and x„+i — +oo). Hence for n > 1 equation (2.3) must be restricted 
in each interval [xk,Xk+i] with the integrals 

Xk + l 

p{x,t)dx, (4.6) 

constrained to a constant value for t > to- This constant is not, in general, equal to one (only the 
sum of these n+1 integrals amounts to one) and, since the separate intervals cannot communicate, 
it will be fixed by the choice of the initial conditions. Therefore, due to the singularities appearing 
in the forward velocity fields «(+)„ for n > 1, we deal with a Fokker-Planck problem with barriers. 
The boundary conditions associated to (2.3) then require the conservation of probability in each 
sub-interval [xk, Xk+i], i.e. the vanishing of the probability current at the end points of the interval: 

[D9.P-^i+)P].k,.k,^=^^ t>to. (4.7) 

To obtain a particular solution one must specify the initial conditions. In particular, we are 
interested in the transition probability density p{x , t\xo , to) , which is singled out by the initial 



condition (3.4), because the asymptotic convergence in L^ of the solutions of equation (2.3) is 
ruled by the asymptotic behaviour of p(a;, t\xQ, Iq) through the Chapman Kolmogorov equation [1] 

/ + 00 
p{x,t\y,to)p{y,t^)dy. (4.8) 

-oo 

It is clear at this point that in every interval [a;fe,a;fc+i] (both finite or infinite) we can write the 
solution of equation (2.3) along the guidelines sketched in Sect. 3. We must only keep in mind 
that in [a;fc,a;fe_|_i] we already know the invariant, time-independent solution 4>n{x) which is never 
zero inside the interval itself, with the exception of the end points Xk and Xk+i- Hence, as we have 
seen in the general case, with the position 

p{x,t) ^ (l3„{x)g{x,t) , (4.9) 

we can reduce (2.3) to the form 

dtg = Cn9, (4.10) 

where >C„ is now the self-adjoint operator defined on [.t^;, a:;/j+i] by 

qn{x)(pix) , (4.11) 



CM^) = ^ 



dip{x) 
r(x)' 



dx 



with 



V? , X (x) v', , s (x) 



r[x)=D>^, ,„(^) = ^^^ + ^Z^. (4.12) 

Equation (4.10) is solved by separating the variables, so that we immediately have 7(t) = e~'^* 
while the spatial part G(x) of g must be solution of 

/:„G(x) + AG(x) == , (4.13) 

with the boundary conditions 

{2DG\x) - «(+)„(x)G(x)] ^^^^^^^ = . (4.14) 

The general behaviour of the solutions of this Sturm-Liouville problem obtained as expansions in 
the system of the eigenfunctions of (4.13) has already been discussed in Section 3. In particular 
we deduce from (3.12) that for the stationary quantum states (more precisely, in every subinterval 
defined by two subsequent nodes) all the solutions of (2.3) always converge in time toward the 
correct quantum solution |0riP- As a further consequence, any quantum solution (^\ defined on 
the entire interval (— oo, +oo) will be stable under deviations from its initial condition. 

5. An explicit example: the harmonic oscillator 

To provide an explicit evolution of the probability and the transition probability densities of 
stochastic mechanics, we consider in detail the example of a harmonic oscillator associated to 
the potential 

n-) = y-V, (5.1) 



with energy eigenvalues 

En^nuj(n + -], n^0,l,2.... (5.2) 

Introducing the notation 



n 

2muj 
the time-independent part of the eigenfunctions (4.1) reads 



'^o = :^ , (5-3) 



M^) - , I- , c-^^/^^° H^ (^) , (5.4) 

where Hn are the Hermite polynomials. The corresponding forward velocity fields for the lowest 
lying levels are: 

W(+)i(x) = 2— --wx, (55) 

with singularities in the zeros of the Hermite polynomials. When n = the equation (2.3) takes 
the form 

dtp ~ uja'^d'^p + ujxdxP + i^p , (5-6) 

and the fundamental solution turns out to be the Ornstein-Uhlenbcck transition probability density 

Po{x, t\xo,to) = —^-j= e-[--"(*)l'/2-'(*) , (i > to) , (5.7) 

o'(i)v27r 

where we have introduced the notation 

a{t) = xoe--^*-*") , a^it) = a^ [l - e-2-(*-*«)] , {t > to) . (5.8) 

The stationary Markov process associated to the transition probability density (5.7) is selected by 
the initial, invariant probability density 

Po(:r:) = ^=e-^/2'^o, (5.9) 



which is also the asymptotic probability density for every other initial condition when the evolution 
is ruled by equation (5.6) (see [1]) so that the invariant distribution plays also the role of the 
limit distribution. Since this invariant probability density also coincides with the quantum one 
0Q ~ iV'oPj the process associated by stochastic mechanics to the ground state of the harmonic 
oscillator is nothing but the stationary Ornstein-Uhlenbcck process. 

For n > I the solutions of (2.3) are determined in the following way. As discussed in the 
previous section, one has to solve the eigenvalue problem (4.13) which can now be written as 

_ |_ G"{x) + (^uj^x^ - Tiiu ^^^) G{x) ^ n\G{x) , (5.10) 

2m \ 2 2 y 
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in every interval [xk, Xk+i] between two subsequent singularities of the forward velocity fields f(+)n. 
The boundary conditions at the end points of these intervals, deduced from (4.14) through (4.5), 
are 

[(/)„G" - (t>nG]^^,^^^^ = . (5.11) 

Reminding that 0„ (but not (f>[^) vanishes in Xk,Xk+i, the conditions to impose are 

Gixk) = Gixk+i) ^ , (5.12) 

where it is understood that for xq and a;„+i we mean, respectively 

lim G{x)^0, lim G{x) ^ . (5.13) 

It is also useful at this point to state the eigenvalue problem in adimensional form by using the 
reduced eigenvalue /i — X/lo, and the adimensional variable x/ao which, by a slight abuse of 
notation, will be still denoted by x. In this way the equation (5.10) with the conditions (5.12) 

becomes 

fx^ 2n + l \ f . ,, 
y (x) — \ n y[x) ~ , 

yi^k) ^y{xk+i) =0. 
If fXm and ym{x) are the eigenvalues and eigenfunctions of (5.14), the general solution of the 
corresponding Fokker-Planck equation (2.3) will be (reverting to dimensional variables) 

X 



p(x, t)^J2 Cme-'^'""Vn(a;)ym — • (5-15) 



The values of the coefficients Cm will be fixed by the initial conditions and by the obvious require- 
ments that p{x, t) must be non negative and normalized during the whole time evolution. Two 
linearly independent solutions of (5.14) are 

y'^' = e-^V%/ (- ^, i; ^) , .(^) = xe-^^/^M (- ^i±^, I ^) , (5.16) 

where M{a, b; z) are the confluent hypergeometric functions. The complete specification of the 
solutions obviusly requires the knowledge of all the eigenvalues /im- 

We consider first the instance n — 1 [xq = — oo, xi = and X2 = +oo), which can be 
completely solved [5]. In this case equation (5.14) must be solved separately for x < and for 
X > with the boundary conditions y(0) = and 

lim y{x) = lim y{x) = 0. (5-17) 

X — * — oo ' X — >+oc ' 

A long calculation [5] shows that the transition probability density is now 

^ g^[a;-Q(t)]^/2(T^(t) _ ^-[x+a(t)]'^/2a'^{t) 

p{x,t\xa,to) = —- _- , (5.18) 

a{t) cr(t)V27r 
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where a{t) and cr^(t) arc defined in (5.8). It must be remarked however that (5.18) must be 
considered as restricted to x > when xq > and to x < when xq < 0, and that only on these 
intervals it is suitably normalized. In order to take into account both these possibilities we can 
introduce the Heavyside function 0(x) so that for every xg 7^ we will have 

Pi{x,t\xo,to) = 0{xxo)—— , , ^— . (5.19) 

^From equation (4.8) we can deduce the evolution of every other initial probability density. In 
particular it can be shown that, with pi(.x) = '/'i(x) 

lim pi(x,t|xo,io) = 2e(xxo) „ ,-— ^ 29(xxo)pi(x) . (5.20) 

Hence, if p{x,tQ) — pin{x) is the initial probability density, we have for t > to 

/+00 
p{x,t\y,to)pin{y)dy 

'Z (5.21) 

Q{xy)ptn{y) dy = r(e; x)pi(x) , 

-00 

where we have defined the function 

r(e;x) = ee(x) + (2-e)e(-x); e^2 pUv) dy ■ (5.22) 

Jo 

When e = 1 (with symmetric initial probability, equally shared on the two real semi-axis) we 
have r(l;x) = 1 and the asymptotical probability density coincides with the quantum stationary 
density pi (x) = (fif (x) . If on the other hand e 7^ 1 the asymptotic probability density has the same 
shape of (j>1{x) but with different weights on the two semiaxes. 

If we consider the higher excited states, the Sturm-Liouville problem (5.14) must be solved 
numerically in each sub-interval. For instance, in the case n = 2 we have xq = —00, xi = —1, 
X2 = 1 and X3 = +00. Considering in particular the sub-interval [—1,1], it can be shown that 
beyond /io = the first few eigenvalues are determined as the first possible values such that 

M(-^.|;i)^o, ,5.3, 

This gives pi ~ 7.44, /xa ~ 37.06, p3 ~ 86.41. 



6. Controlled evolutions 

In this Section we move on to implement the program outlined in the introduction, that is to exploit 
the transition probabilities of Nelson stochastic mechanics to model controlled quantum evolutions 
to arbitrarily assigned final quantum states. We start by observing that to every solution p{x, t) of 
the Fokker-Planck equation (3.1), with a given W(+)(x,i) and constant diffusion coefficient (2.11), 
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we can always associate the wave function of a quantum system. To this end, it is sufficient to 
introduce a suitable time-dependent potential. 

Let us take a solution p(x, t) of the Fokker-Planck equation (3.1), with a given vi^\{x, t) and 
a constant diffusion coefficient D: introduce the functions R{x, t) and W{x, t) from 

p{x,t) = R'{x,t), V(^+){x,t)^d^W{x,t), (6.1) 

and remind from (2.12) that the relation 

mv(+)^d,S + n^^d,S+\^^d,[s+^\np^ (6.2) 

must hold, where p is an adimensional function (argument of a logarithm) obtained from the 
probability density p by means of a suitable and arbitrary dimensional multiplicative constant. 
If we now impose that the function S{x, t) must be the phase of a wave function as in (2.8), we 
immediately obtain from (6.1) and (6.2) the equation 

S{x,t) ='mW{x,t) - -\np{x,t) - e(t) , (6.3) 

which allows to determine S from p and U(_|_) (namely W) up to an additive arbitrary function of 
time 0(t). However, in order that the wave function (2.8) with R and S given above be a solution 
of a Schrodinger problem in quantum mechanics, we must also make sure that the Hamilton- 
Jacobi-Madelung equation (2.9) is satisfied. Since S and R are now fixed, equation (2.9) must be 
considered as a relation (constraint) defining the controlling potential Vc, which, after straightfor- 
ward calculations, yields 

fi2 t ravj.s 
V,{x,t)^—dlhip+-{dt\np + Vi^+)dA^p) ^-mdtW + e. (6.4) 

Of course if we start with a quantum wave function tp{x^ t) associated to a given time-independent 
potential V{x) and we select as a solution of (2.3) exactly p — ji/ip, then the formula (6.4) always 
yields back the given potential, as it should. This can be explicitly seen (to become familiar 
with this kind of approach) in the examples of the ground state and the first excited state of the 
harmonic oscillator potential (5.1), by choosing respectively in equation (6.4) 6{t) ~ hujt/2 and 
6{t) = 3hLut/2, which amounts to suitably fix the zero of the potential energy. 

On the other hand the nonstationary fundamental solution (5.7) associated to the velocity field 
f(+)o(a;) of (5.5) for the case n = (we put fo = to simplify the notation) does not correspond to 
a quantum wave function of the harmonic oscillator whatsoever. However it is easy to show that, 
by choosing 

and the time-dependent controlling potential 

2 _9 ___ .9_9 ___ T^T 



v;(.,t).^ 



X — Oi\t) 

^(iT 
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(t -^ +oo) , (6.6) 



with a{t), a{t) and (Tq defined in equations (5.8) and (5.3), we can define a quantum state, i.e. 
a wave function -tpcix^t) solution of tlie Sclirodinger equation in the potential (6.6). At the same 
time Tpc is associated to the transition probability density of the form (5.7) which is its modulus 
squared. Of course the fact that for t -^ +oo we recover the harmonic potential is connected to 
the already remarked fact that the usual quantum probability density is also the limit distribution 
for every initial condition and in particular also for the evolution (5.7). In the case n = 1, with 
the W(_|_)i(a;) as given by equation(5.5) and the transition probability density (5.19), we define 



Tix)^---, (6.7) 

tanha; 



and then we choose 



m--[^^-^^-l)^-,n^, (t-+oo), (6.8) 



so that we have the following time-dependent controlling potential (for every x =/= 0): 



4 



(t^ + 



^_ Og^/a^ 



Amx^ 



,' xa 



(6.9) 



2 ' 
The limit t -^ +oo must be obviously intended in a physical sense, i.e. for times much longer 

than X^ , the largest characteristic time of decay in the expansion (3.12). In this particular case 

Ai = uj. In fact here too the asymptotic potential is the usual one of the harmonic oscillator, but 

it must be considered separately on the positive and on the negative x semiaxis, since in the point 

x — Q a singular behaviour would show up when i ^ 0. This means that, also if asymptotically 

we recover the right potential, it will be associated with a new boundary condition in a; = since 

the system must be confined on one of the two semiaxes. 



7. Modelling transitions 

Given any couple [p, «(+)) associated to a Fokker Planck equation, the possibility of promoting it 
to a solution of a Schrodinger problem by a suitable controlling potential Vc{x, t) enables to model 
quantum evolutions driving, for instance, the probability density of a given quantum stationary 
state to another (decays and excitations). Moreover, an immediate generalization of this scheme 
might open the way to modelling evolutions from a given, arbitrary quantum state to an eigen- 
function of a given observable. Besides other applications, this is something which could become 
a starting point for building simple models of the measurement process, where one tries to dy- 
namically describe the wave packet collapse [10]. As a first example let us consider the transition 
between the invariant probability densities associated to the ground and the first excited state of 
the harmonic oscillator potential (5.1): 

""f^ (7.1) 

Pl(x)=0?(x) = ^^e--V2.o^ 



^0' 
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If we choose to describe the decay (jji -^ (j}f) we may exploit the Chapman-Kolmogorov equation 
(4.8) with the transition probabihty density (5.7), and the initial probability density pi{x). An 
elementary integration shows in this case that the resulting evolution takes the form (tg = 0) 

pi^o(a:,t)=/3'(t)po(a;)+7'WPiW, (7.2) 

where 

/?2 (t) = 1 - e-2-* , ^{t) = e-"* . (7.3) 

Recalling V(^-^q[x) as given in (5.5) and the evolving probability density (7.2), and inserting them 
in equation (6.4) we obtain the following form of the controlling potential: 

Tnijj T 

Vc{x, t) = 2huuU{xlao- /3/7) , (7.4) 



where 



The parameter 



^(-^^)^ (52 +.2). ■ (7-5) 






^W = ^T7T=c 1 (7.6) 



is such that 62(o+) = and b'^{+oo) — +oo. Thus U goes to zero as i ^ +oo for any x, and as 
t ^ 0+ is 1 for every x, except for a negative singularity in a; = 0. As a consequence, while for 
t -^ +00 the controlling potential (7.4) simply tends to the potential (5.1), for i ^ 0+ it presents 
an unessential shift of —2huj in the zcroth level, and a deep negative singularity in x = 0. 

The singular behaviour of the controlling potential at the initial time of the evolution is a 
problem connected to the proper definition of the phase function S. In fact, from (6.3) we have: 



S{x,t) = "2 ^^ 



/32(i) + ^7'(t) 



T^' (") 



Six,o+) = ---\n::^. (7.. 



so that in particular we have 

h 

In- 

We would instead have expected that initially the phase be independent of x as for every stationary 
wave function. This means that in the constructed evolution S{x, t) presents a discontinuous 
behaviour for t ^ 0^. The problem arises here from the fact that we initially have a stationary 
state characterized by a probability density pi{x) and a velocity field v/^\i(x), and then suddenly, 
in order to activate the decay, we impose to the same pi to be embedded in a different velocity 
field v/^\q(x) which drags it toward a new stationary po{x). This discontinuous change from W(+)i 
to f (+)o is of course responsible for the remarked discontinuous change in the phase of the wave 
function. We have therefore modelled a transition which starts with a sudden, discontinuous kick. 
To construct a transition that evolves smoothly also for t ^ 0+ we should take into account a 
continuous and smooth modification of the initial velocity field into the final one. This requirement 
compels us to consider a new class of Fokker-Planck equations with time-dependent forward 
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velocity fields U(_|_)(x,t). In particular, to achieve the proposed smooth controlled decay between 
two stationary states of the harmonic oscillator, we should solve an evolution equation with a 
continuous velocity field v/^\{x,t) which evolves smoothly from vt_^_\l(x) to vt_^_\Q(x). Clearly, the 
smoothing procedure can be realized in several different ways and the selection must be dictated 
by the actual physical requirements and outputs one is interested in. A suitable smoothing for our 
transitions which leads to manageable equations still has to be found; however in the following 
Section we will study a problem in which the smoothness of the evolution is a priori granted. 



8. Smooth transitions: coherent and squeezed wave packets 

As anticipated at the end of the previous Section we will now consider an instance of controlled 
evolution that does not require an extra smoothing procedure for the driving velocity field, i.e. 
the transition between pairs of coherent wave packets. In particular we will consider both the 
transition from a coherent oscillating packet (coherent state) to the ground state of the same 
harmonic oscillator, and a dynamical procedure of squeezing a coherent wave packet. 

To this end we will recall a simple result [1] which indicates how to find the solutions of a 
particular class of evolution equations (2.3) which includes the situation of our proposed examples. 
If the velocity field of the evolution equation (2.3) has the linear form 

V(+){x,t)=A{t) + B{t)x (8.1) 

with A{t) and B{t) continuous functions of time, then there are always solutions of the Fokker- 
Planck equation in the normal form Af[fi{t), i^(i)), where ^{t) and i'{t) are solutions of the differ- 
ential equations 

/i(t) - B{t)ii{t) = A{t) ; i>(i) - 2B{t)v{t) = 2D (8.2) 

with suitable initial conditions. The first case that we consider is the coherent wave packet with a 
certain initial displacement a: 



^(^' ^) = ( ^:^ ) ^^P 



1 \ r {x — a cos bjt)"^ / Aaxs\nu)t — a? sv[\2Lot utt 

i^g \ 8^^ +Y 



^0 



(8.3) 



whose forward velocity field reads 

U(+) (x, t) — aLu{cosLut — sinLut) — lux . (8-4) 

The field (8.4) is of the required form (8.1) with A{t) = aLu{cosLLjt — sinujt) and B{t) ~ ~uj, while 
the configurational probability density is 

p{x,t) — |V'(a;,t)p = po{x — acosujt) , (8.5) 

with po is the one of the ground state of the harmonic oscillator given by (7.1). It is easy to show 
that when B{t) = —lu, as in the case of the wave packets we are considering, there are coherent 
solutions of (2.3) with u{t) = ctq of the form A/'(/i(t), 0-0)1 i-C- of the form 

p{x,t)=po{x-fi{t)). (8.6) 
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Now the time evolution of such coherent sohitions can be determined in one step, without imple- 
menting the two-step procedure of first calculating the transition probability density and then, 
through the Chapman-Kolmogorov equation, the evolution of an arbitrary initial probability den- 
sity. On the other hand if we compare (5.5) and (8.4) we see that the difference between U(+)o 
and W(+) consists in the first, time-dependent term of the latter; hence it is natural to consider the 
problem of solving the evolution equation (2.3) with a velocity field of the type 

V(,) [x, t) = A{t) — ujx , 

(8.7) 
A(t) = aLjj(cosLjjt — shiLjjt)F{t) , 

where F{t) is an arbitrary function varying smoothly between 1 and 0, or viceversa. In this case 
the evolution equation (2.3) still has coherent solutions of the form (8.6) with a /i(t) dependent on 
our choice of F{t) through equation (8.2). 

A completely smooth transition from the coherent, oscillating wave function (8. .3) to the 
ground state (j)^ (5.4) of the harmonic oscillator can now be achieved for example by means of the 
following choice of the function F{t): 

N 



F(t) = l-(l-e-^*)^ = ^(-l)'=+i(^) 



fc=i 

where 

IniV 
n^ , LOk = kQ.- T>0, N>2. (8.9) 

T 

In fact, a function F{t) of this form goes monotonically from F{0) = 1 to F{+oo) = with a flex 
point in r (which can be considered as the arbitrary instant of the transition) where its derivative 
F'{t) is negative and grows, in absolute value, logarithmically with N . The condition that the 
exponent N >2 also guarantees that F'{0) — 0, and hence that the controlling potential Vc{x,t) 
given in equation (6.4) will continuously start at i = from the harmonic oscillator potential (5.1), 
and asymptotically come back to it for t -^ +oo. Finally, the phase function S{x, t) will too change 
continuously from that of ip given in (8.3) to that of the harmonic oscillator ground state V'o- A 
long calculation yields the explicit form of the controlling potential: 

V,{x,t)^mLO^—-mujaxY,{-lf+Hmk{t)oJkC-'^'''~WuuJe-'^'] , (8.10) 

fc=i ^ ^ 

where 

, , . 21^^ &\nu)t — wt cos tot 

Uk[t) — smujt + 



[cjk — Uj)'^ + Cij2 



(8J1) 



The parameters r and N, apart the constraints (8.9), are free and can be fixed by the particular 
form of the transition that we want to implement, according to what specific physical situations we 
are interested in. Finally we remark that, in a harmonic oscillator, the transition between a coher- 
ent, oscillating wave packet and the ground state is a transition between a (Poisson) superposition 

17 



of all the energy eigenstates to just one energy eigenstate: an outcome whicli is similar to that of an 
energy measurement, but for the important difference that here the result (the energy eigenstate) 
is deterministically controlled by a time-dependent potential. The controlled transition that we 
have constructed does not produce mixtures, but pure states (eigenstates) and may be considered 
a dynamical model for one of the branches of a measurement leading to a selected eigenvalue and 
eigenstate. 

Until now we have considered transitions between gaussian wave packets with constant width. 
However it is also of great interest to discuss the case of controlling potentials able to produce a 
wave-packet evolution with variable width: a kind of controlled squeezing of the wave-packet. This 
could be very useful in such instances as the shaping of the Gaussian output in the manufacturing 
of molecular reactions, or in the design of focusing devices for beams in particle accelerators, in 
which the width of the bunch must be properly squeezed. We will discuss now a simple case which 
shows also that, in the particular conditions chosen, it is also possible to avoid the integration of 
the differential equations (8.2). 

Let us remember that when the forward velocity field has the form (8.1) the Fokker-Planck 
equation (2.3) always possesses Gaussian solutions of the form 

^^-[x-^,{t)f/2,y(t) 

if /i(t) and ^{t) are solutions of (8.2). We plan now to describe evolutions of the quantum state 
(2.8) such that I) both Vc{x,t) and S(x,t) be continuous and regular at every instant, and II) the 
variance i'{t) satisfy the relations 

i/(-oo) = cr^, z/(+oo) = cr^. (8.13) 

In practice this means that, if for example we require for the sake of simplicity fi(t) = at every 
time, we will describe a transition from the ground state of an harmonic oscillator with frequency 
loq =^ D/a^ to the ground state of another harmonic oscillator with frequency uji = D/af. It 
is convenient to remark here that this very simple transition cannot be achieved by means of an 
arbitrary time-dependent potential Vc{x,t), given that it goes from mio^x'^ /2 for t -^ — oo to 
mojix'^ /2 for t -^ +oo. The intermediate evolution, indeed, when not suitably designed, would 
introduce components of every other energy eigenstate of the final harmonic oscillator which will 
not, in general, asymptotically disappear. 

Let us recall here that the relevant quantities are the phase function 

S{x,t) ^mW{x,t) -mD\np{x,t) - 9{t) (8.14) 

(where 9{t) is arbitrary and, from (6.1) and (8.1), W{x,t) ~ A{t)x + B{t)x'^ /2), and the controlling 
potential 

Vcix,t) =mD^dl\iip + mD{dt\iip + v^+)d^lnj5) ^ - mdtW + 6 . (8.15) 

18 



Both these two functions are determined from the knowledge of forward velocity field W(_|_) (x, t) 
and of the adimensional density p{x,t) = (7Qp(x,t). However in this coherent evolution it will not 
be necessary to integrate the differential equations (8.2) to obtain an explicit form of S and Vc- 
Indeed, since A{t), B{t) and \np{x,t) can be expressed through (8.2) in terms oi iJ,{t), i'{t), fi{t), 
t'(t) and D, it is a straightforward matter to show that the phase is of the general form 

S{x, t) = y [n{t)x^ - 2U{t)x + A{t)] , (8.16) 

with 



2v{t) 



^ * = D —7— +Dln 2 , 

h'[t) ctq m 



while the controlling potential reads 



V^{x, t) = y [uj^{t)x^ - 2a{t)x + c{t)] , (8.18) 



where 

AD^ - 2iy{t)i>{t) + i^^t) 



«W = JTJTT^ ' (8-19) 



C{t) 



AD'^p{t) + Av'^{t)il{t) - 2p{t)v{t)v(t) + p(t)P(t) 

8D^p^(t) - ADv{t)v{t) - 8D^v{t) - {2v{t)p{t) - p{t)v{t) + 2Dp{t)) 20{t) 



4i.2(i) 



We can simplify our notation by imposing that fi{t) = (and hence /i(t) = /i(t) = 0) for every t, 
obtaining 

S{x,t)^'^[n{t)x^ + A{t)], (8.20) 



with 



2v{t) ' 

A(t)^i.ln^-« 
Cq m 

The controlling potential too reduces to 



(8.21) 



Vc{x, t) = y [u\t)x' + c(i)] , (8.22) 



where 



2 _ AD^-2u{t)i){t) + i^Ht) 
'"(*^- Ai^) ' 

2e{t) iA^)m + 2lM^ 
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.23) 



Hence the evolution is completely defined, through the four functions ri(t), A(t), uj'^{t) and c{t), 
by 9(t) and i^(t). It is expedient in particular to choose 

.(t)^ — ln-^ + -^. (8.24) 

In this way 

A{t)^-^, (8.25) 

because from (8.13) we have i>(±oo) = so that (see (8.20)): 

b[X,t) ^ 2 : t ^ — OO , 

1, (8.26) 

S[X,t) ^ 2 ' t ^ +00 . 

^\ 

This was to be expected from the fact that mD^ ja^ = huJo/2 and mD^ /a^ = huJi/2 are the energy 
eigenvalues of the ground states of the two harmonic oscillators. Moreover from the choice (8.24) 
also follows that c^{ztoo) = so that the controlling potential (8.22) will show no asymptotical 
extra terms with respect to the initial and final harmonic potentials. 

In order to completely specify the controlled evolution we are now left with the determination 
of the form of i'{t). li b = cti/cto, then take 

-W = -0^ (yT^TTT^) . (->0), (8.27) 

so that the transition happens around the instant t = and t controls its velocity. We thus obtain 
the explicit expressions for the four functions (8.21) and (8.23): 

n{t) 

A(t) 

^'-^^^^ a2 (6 + e-*/-)3 , 

Their form is displayed in the Figures 1-4, where to fix an example we have chosen the values 
T=l,& = 2, a;o = l and cro = 1- As it can be seen in this case the behaviour of the potential 
time-dependent parameters is not trivial even for the very simple squeezing of a static gaussian 
wave packet from a given variance to another. How to precisely follow this time dependence in a 
stable way will be the argument of a forthcoming paper, as discussed in the next section. 



b- 1 e 


-t/T 








r (6 + e-*A 


■)(l + e-*/- 


0' 






2DH /1 + e- 


-"f 








<jf) \b + c- 


-„r) • 








D^ /1 + e-*/^ 

(T^ \b + C-t/^ 
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(1 + c 


t/r(l_c-*/^) 
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(8.28) 
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9. Conclusions and outlook 

We have shown how to treat the typical inverse problem in quantum control, i.e. that of 
determining a controlling potential for a given quantum evolution, in the framework of Nelson 
stochastic mechanics. In this way we have been able to determine the general characteristics of 
controlled evolutions between assigned initial and final quantum states. The techniques of solution 
and the relation between the transition probabilities, phase functions and controlling potentials 
have been discussed on general grounds. Detailed, explicit calculations have also been shown in 
the paradigmatic test arena provided by the harmonic oscillator. 

Further extensions of the method outlined in the present paper are currently under study. One 
immediate application to be faced is the generalization of the analysis performed for the harmonic 
oscillator to anharmonic systems. The difficulty to be faced on the way toward this aim is that 
one is in general forced to deal with approximate quantum wave functions, as in the case of the 
quartic oscillator. Therefore the controlled evolution must be supplemented by a suitable feedback 
mechanism ensuring that the error initially made in choosing a certain initial approximate state 
does not grow during the controlled time evolution. One extremely interesting application would be 
the description of a controlled evolution driving initial approximate quantum states of anharmonic 
systems to stable wave packets generalizations of the coherent states of the harmonic oscillator 
[15]. Besides the obvious interest in several areas of quantum phenomenology, the above is also of 
great potential interest in discussing the control and the reduction of aberrations in quantum-like 
systems, i.e. deviations from the harmonic evolution that are detected in systems like charged 
beams in particle accelerators. 

Another very interesting future line of research that has been left virtually unexplored in 
the present paper is the introduction of optimization procedures. We have barely touched upon 
this problem when discussing the smoothing of the controlled transitions. Optimization of suitable 
functionals, chosen according to the kind of physical evolution one needs or desires to manufacture, 
would provide a powerful criterion of selection among the different possible smoothed evolutions. 
Instances of functionals to be optimized during the controlled dynamics that come naturally to 
mind are the uncertainty products of conjugate observables (to be optimized to a relative minimum 
under the constraint of Schrodinger dynamics [16]), or the relative entropy between the initial and 
final states. But many more can be imagined and devised, according to the nature of the physical 
problem considered. 

One last, but important consideration is in order. When we implement a controlled evolution 
by means of a suitable controlling potential we must also bear in mind that in practice small 
deviations away from the designed potential and from the desired wave function are always possible. 
In general such deviations are not subsequently reabsorbed but rather tend to drag the state away 
from the required evolution. Hence to really control these quantum evolutions it will be very 
important to study their stability under small deviations and perturbations: this is of crucial 
importance from the standpoint of confronting the formal, theoretical scheme with the practical 
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applications. Work is currently in progress in all the above mentioned extensions of the present 
research, and we plan to soon report on it. 
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FIGURE CAPTIONS 

Figure 1. The parameter 51(t) in the phase function S{x,t). 
Figure 2. The parameter A(i) in the phase function S{x,t). 
Figure 3. The parameter w^(i) in the potential Vc{x,t). 
Figure 4. The parameter c{t) in the potential Vc{x,t). 
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